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We computed the bottom quark forward-backward asymmetry at the Tevatron in the Standard 
Model and for several new physics scenarios. Near the Z-po\e, the SM bottom asymmetry is domi- 
nated by tree level exchanges of electroweak gauge bosons. While above the Z-po\e, next-to-leading 
order QCD dominates the SM asymmetry as was the case with the top quark forward-backward 
asymmetry. Light new physics, Mjvp ^ 150 GeV, can cause significant deviations from the SM 
prediction for the bottom asymmetry. The bottom asymmetry can be used to distinguish between 
competing NP explanations of the top asymmetry based on how the NP interferes with s-channel 
gluon and Z exchange. 
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INTRODUCTION 

Measurements [1-3] of the forward-backward 
asymmetry in top quark pair production (A^pg) by 
the CDF and D0 collaborations at the Tevatron 
have attracted a lot of attention recently. At high 
invariant mass, the CDF measurement ^"^(-^^tt — 
450 GeV) = 0.295 ± 0.058(stat.) ± 0.031(syst.) is 
approximately Bcr away from the Standard Model 
(SM) prediction, 0.100 ± 0.030 [ ]. In addition, 
CDF observes that has an approximately lin- 
ear dependence on both the invariant mass and the 
magnitude of the rapidity difference (|Aytf|) of the 
tt pair with slopes that are more than 2a away 
from the SM prediction. 

Soon after CDF reported evidence for a mass- 
dependent tt asymmetry, it was realized [4-6] that 
measuring the forward-backward asymmetry in 
bottom quark production (A^pg) may provide in- 
sight into the source of the tt asymmetry. Any new 
physics (NP) explanation of Ap^ that respects 
SU{2)]^ or custodial symmetry will in general also 
create an asymmetry in bb production. The CDF 
collaboration is in the process of measuring the bb 
forward-backward asymmetry, and has stated [ , ] 
how it is binning the data and how sensitive it ex- 
pects to be to a potential signal. However, A'^g 
will likely be more difficult to measure than A"^. 
Among the reasons for this are that gluon fusion, 
which does not produce an asymmetry, is responsi- 
ble for >, 90% of bottom quark production at the 
Tevatron. In addition, the bb asymmetry is mea- 
sured by selecting dijet events containing a soft 
muon, and relating the charge of the muon to the 
charge of the b that produced it [ ]. However, as 
explained in more detail in [5, (>], B—B mixing and 
cascade decays will partially wash out the correla- 
tion between the charge of what is detected and 
the charge of the bottom quark that produced it, 
and this needs to be accounted for. [03] 



In this Letter, we computed the bottom quark 
forward-backward asymmetry at the Tevatron in 
the SM and for several NP scenarios. It is nec- 
essary to know the SM prediction in order to de- 
termine whether or not any NP can possibly be 
present. Since a small asymmetry is expected in 
the SM, ApB provides an excellent window to ob- 
serve NP. An interesting difference between the 
bottom and top quark asymmetries is that the Z- 
pole is in the signal region for the bb asymmetry. 
This leads to tree level exchanges of electroweak 
gauge bosons dominating the SM contribution to 
ApB near the Z-pole, as well as the opportunity for 
there to be significant interference effects between 
NP and tree level Z exchange. 



STANDARD MODEL CALCULATION 

The definition of the forward-backward asym- 
metry in heavy quark production we use is 
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Here Ay is the difference in the rapidity of the 
quark and anti-quark. Ay = yq — yQ, and is invari- 
ant under boosts along the collision axis. A frame 
dependent asymmetry may also be defined using 
yq instead of Ay as the discriminating observable. 
Leading order (LO) QCD is completely symmet- 
ric with respect to Ay, and thus does not generate 
an asymmetry. Starting with next-to-leading or- 
der (NLO) QCD, contributions to the asymmetry 
as an expansion in powers of can be written 
schematically as 
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Analytic formulae for the 0{as) and 0{a) terms 
of ApB are given in ['', ']. These results are 
based on analogous calculations [10, 11] for the 
e~e'^ — >■ 7* — >■ fJ-~ fJ,'^ asymmetry. Prior results 
on the QCD asymmetry also exist [12-14]. The 
0{a^/a^) term for was computed in [1-j]. 

Electroweak Sudakov corrections are shown in [Ifi] 
to increase the 0{as) contribution to the inclu- 
sive by a factor of 1.07. While the iVi and 
Di terms in (2) are known completely and have 
been studied [17-26] in depth, N2 is only partially 
known [27, 28]. [64] Since it would be inconsistent 
to include the NiDi/Dq term in our calculation 
without the N2 term, we drop the 0{a^) contribu- 
tion to ApB- To account for this neglect of higher 
order terms, we assign an uncertainty to our cal- 
culation of 30% of the 0{as) contribution, origi- 
nating from oigDi 0.3I?o- 

Our calculation was done by convolving the ana- 
lytic formulae of [9, 15] with MSTW2008 PDFs [29] 
using the deterministic numeric integration algo- 
rithm Cuhre from the CUBA library [M)]. LO 
PDFs were used for all LO cross sections, and sim- 
ilarly NLO PDFs were used for NLO cross sec- 
tions, as is set by the MSTW2008 best-fit value, 
as{Mz) = 0.139, 0.120, for the LO and NLO PDFs 
respectively. We fixed ^lji — = Mz and 
nif =4. The other numeric values employed in 
this analysis were: mf, ~ 4.7 GeV, Mz — 91.1876 
GeV, Vz = 2.4952 GeV, a{Mz) = 1/128.93, and 
sm^Ow = 0.231. 

To mimic CDF's analysis [7] we required the 
bh pair in our calculation to have a maximum 
acoUinearity of 5 = tt — 2.8 radians and each h 
to have a minimum energy of 15 GeV. The phase 
space that is available to the gluon in the hhg fi- 
nal state is discussed in [■>!]. We used the soft 
gluon approximation for the amplitude to hhg; the 
error so induced, estimated in ['! '], is smaller than 
other errors discussed below. An additional cut, 
1 2/6 6 1 ^ Ij was made as this is roughly the max- 
imum rapidity at which CDF can detect muons 
and 6-tagging quickly loses efficiency outside of this 
region [-!, •>]. We found the 0{a) corrections de- 
crease the contribution of ©(as) to A^p^ by 3-11%, 
depending on the bin. However, we neglect this 
0{a) contribution as it is mostly canceled by the 
increase in A^p^ due to electroweak Sudakov ef- 
fects [10], and the sum of the two effects is small 
compared to the uncertainty in the total contribu- 
tion. The flavor excitation process, qg — qhh, as 
well as ^-channel W exchange were also neglected 
as they are numerically small [9, 15]. 

Our results for the 0{a^ /a1) and 0{as) contri- 
butions to binned A^pp, are shown in table L In the 



second and third columns the uncertainty is due 
to varying /i^ — fJ-F from A/2/2 to 2Mz. In the 
fourth column the first uncertainty is due to ne- 
glect of higher-order terms, and the second is the 
combined scale uncertainty. The uncertainty in the 
0(0^/0^) contribution to A^pg is larger than the 
0{ag) term because the extra power of ag makes 
it more sensitive to the choice of scales and PDFs. 

Based on CDF's expected sensitivities [ ] and 
assuming the Standard Model (and the measure- 
ments follow a Gaussian distribution), CDF should 
be able to exclude A%{75 < A4JGeV < 95) = 
at the 2.5cr confidence level (CL). Although the 
central value for the asymmetry in the > 130 
GeV invariant mass bin is slightly larger than the 
75 - 95 GeV bin, CDF should only be able to ex- 
clude Afp,{130 < M^-JGeY) = at the 1.3cr CL. 
The likelihood of excluding zero asymmetry in the 
95 — 130 GeV invariant mass bin is comparable 
to the likelihood in the > 130 GeV bin. In the 
SM, all the other (mass or rapidity) bins should 
be consistent with zero at the Icr level based on 
experimental uncertainties. 

LO event generators can predict the 0{a'^ /a^) 
contribution to the asymmetry. MadGraph 
5.1.5.5 [•■' ] with CTEQ6L1 PDFs [ '-S] gives 
Afp,{75 < MfcJGeV < 95) = (2.26 ± 
0.32(stat.)to;74(scale))%, in good agreement with 
our calculation. 

It has been suggested [5, 6] that measuring the 
charm quark forward-backward asymmetry at the 
Tevatron (A'^p^) and the bottom quark charge 
asymmetry at the LHC (Aq) may also provide in- 
sight into the origin of the Apg anomaly. We com- 
puted SM asymmetries of a few percent in suitably 
chosen kinematic regions for both A'^g and A^. 
While the central values for these asymmetries are 
comparable to those of A^^, it is unlikely that 
these asymmetries will be observed any time soon 
in the absence of NP. For A'^g, c-tagging is less 
efficient than 6-tagging. For A'^, the kinematic re- 
gions where the asymmetry becomes a few percent 
have small production cross sections, and will re- 
quire the LHC to run for at least a year at 14 TeV 
to collect enough data for the SM asymmetry to 
be statistically distinguishable from zero. 

NEW PHYSICS SCENARIOS 

Many new physics models have been 
proposed [34-43] as explanations of the 
anomalously['i^] large tt forward-backward 
asymmetry. For the stringent constraints that 
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Bin 0(QVag) 0(q,) A%[%] 

35 < M^s/GeV < 75 
75 < Mjs/GeV < 95 
95 < Mjj/GeV < 130 
130 < M^-JGeY 

0.0 < \ Ay^i\ < 0.5 Omtllfi 0.166tg:;];5i 0.22 ± 0.05 Ig;! 

0.5 < I Ay,5| < 1.0 0.111^2 0A02t°o°w 0-52 ± 0.12 1^;! 

1.0 < |Ay,-,| < 2.0 Omtrol 0.444lg °lg 0.53 ± 0.13 Ig;! 

Inclusive 0-09to:o2 0-389tg;go9 0.48 ± 0.12 tp.! 



9 cc+O.SO 
Z.OJ_(j 73 



0.15 



[-0.05 
-0.04 



n 00(5+0.007 

q64+° °°5 

1 c,qts+oooi 
i.oyo_o.oo2 

O.U^O_Q 032 



0.34 ± 0.10 
3.51 ± 0.29 tg-^g 
2.22 ±0.48 1°;?9 
3.77 ± 1.09 tS S? 



TABLE I: The ©(a^/a^) and 0{as) contributions to Afg in various bins. 



these models must overcome see [44--'')l]. Pre- 
dictions for A''pg in the context of various NP 
scenarios have already been made in [G, 43, 44, o2]. 
We expanded on these works by taking into ac- 
count the resonance effects of the Z, and limiting 
ourselves to the energy regime accessible at the 
Tevatron. In particular, we are interested in 
seeing if the NP contribution to can be large 
enough to be distinguishable from the SM predic- 
tions we computed above based on the expected 
sensitivities given in [7]. Any NP in the bottom 
sector must not spoil the agreement between the 
SM and precise measurements of flavor changing 
decays and meson mixing observables such as 
Br(6 — >■ S-I-7) and B — B mixing. These and other 
constraints, such as same-sign top production, are 
more easily satisfied in flavor symmetric models 
in which the NP particles form complete repre- 
sentations of the quark global flavor symmetry 
group, Gf = SU{3)u^ x 5i7(3)i5^ x SU{3)q,. 
Furthermore, the flavor symmetry guarantees a 
definite relationship between Ap^ and A^p^. We 
consider three different models, a light, broad 
axigluon (G"), a scalar weak doublet {((>), and 
an SU{3)q^ octet of electroweak triplet (EWT) 
vectors (V); see table II. 

It is convenient to split the contributions to the 
forward-backward asymmetry into two terms 

Afb=A'pb+A^^b- (3) 

A^pg contains the 0{as) contribution to A^pg, and 
can be obtained from table I. The 0{a) contribu- 
tion to the asymmetry could also be included in 
A\pg, but we neglect it in what follows. On the 
other hand, A^pg contains the SM 0{a^/a^) con- 
tribution to the asymmetry as well as contributions 
from NP. This includes both pure NP contribu- 
tions and interference between NP and tree level s- 
channel gluon and Z exchange. We calculated A^pg 
using FeynRules 1.6.1 [53] to implement the NP 
models in MadGraph 5.1.5.5 [■^>"] including elec- 



troweak processes (QED=2). For Afg, 10^ events 
were generated for a given set of parameters using 
the CTEQ6L1 [ ] PDFs with the renormalization 
and factorization scales set to m*. For A''pg, 10^ 
events were generated for each mass bin for a given 
set of parameters with = jip = Mz- As was 
the case for the SM analysis, a cut was placed on 
the rapidity of the bottom quarks, \yi, i\ < 1. 

Predictions for the binned tt and bb asymmetries 
from the NP models are shown in the left and right 
columns of figure 1 respectively. Overflow is in- 
cluded in the rightmost bins. The widths of the 
axigluon and the EWT vectors were chosen to be 
10% of their masses. For the scalars, the natural 
width to quarks was used. Axigluon benchmark 
points were taken from table I of [4"i]. Benchmark 
points for the (j) and V models were chosen based 
on adding approximately 10% to the inclusive tt 
asymmetry, having a roughly linear dependence of 
Apg on Mj(, and adding (or subtracting) less than 
1 pb from the ti production cross section at the 
Tevatron. 

We have given three classes of models that can 
accommodate Apg and produce a Apg that is dis- 
tinguishable from the SM prediction. However, 
this is not generally the case. For example, a fla- 
vor octet, EW singlet model (Vj' (T^)^' ^ Vf^Sf 
in table II), can accommodate A*pg without caus- 
ing any significant deviations from the SM predic- 
tions because it only produces bb from dd initial 
states whereas the other models involve uu initial 
states. All three models considered can interfere 
with gluon exchange, (j) and V can also interfere 
with the Z, which dominates the NP contribution 
to A^pg in the Z-pole bin for these models. 

In addition to the Apg anomaly, there is the 
longstanding puzzle of the bb forward-backward 
asymmetry at LEPl, Apg\ which is 2Aa below 
the SM value [ ]. Furthermore, the ratio of the 
partial width Z ^ bb to the inclusive hadronic 
width, Rt, is 2.3a above the SM prediction [55]. 
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Case SM Gf Relevant Interaction Ref. 



G' (8,l)o 


(1,1,1) 


ga {UR(f' Ur + DR(f' Dr - Ql(J!' Ql) 


[36, 37] 


^ (1, 2)1/2 


(3,1,3) 


X [(fP II Vtb Ur + (j)" II Ur) + h.C. 


[38] 


V (1,3)0 


(1,1,8) 




[39, 40] 



TABLE II: The gauge and flavor representations for the models under consideration. Tq and are 

generators of SU{3)q^ and SU{2)l respectively. 
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: Predictions for the binned A^pg (left) and A^pg (right) from the axigluon (top), scalar weak 
(middle), and flavor octet vector (bottom) models. SM predictions are in orange. In black are 
CDF's measurements [ ] and expected sensitivities [ ] for Apg and A'^pg respectively. 
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Assuming only the bottom quark's coupling to 
the Z is modified, the value of bgfih which pro- 
vides the best-fit to the EWPD collected at LEP 
is 0.016 ['(i], which is more than 20% of the LO 
SM coupling. See [41, 42] for attempts to simulta- 
neously explain J^p^ and ^p'g ■ We verified that 
loop corrections induced by the and V models 
considered above, and the tree-level V — Z mixing 
of [""], cannot produce this large of an effect while 
being consistent with constraints from atomic par- 
ity violation [48]. Prospects for measuring 66 and 
ti asymmetries at future linear colliders are exam- 
ined in ['"1. 



CONCLUSIONS 

In summary, we computed J^pg in the SM and 
for several NP scenarios, carefully accounting for 
the Z-pole, which is in the signal region for the 66 
asymmetry. The largest SM contribution to ^pg 
near the Z-pole comes from tree level exchanges of 
Z and 7*. While at higher invariant mass, NLO 
QCD dominates the SM asymmetry. Light NP, 
-Mnp ~ 150 GeV, is needed to generate a 66 asym- 
metry, which CDF would be able to distinguish 
from the SM. J^p^ can be used to distinguish be- 
tween competing NP explanations of J^p^ based on 
how the NP interferes with s-channel gluon and Z 
exchange. 

We thank Aneesh Manohar for advice regarding 
the SM calculation, Jesse Thaler for an enlighten- 
ing discussion, and Michael Trott for useful discus- 
sions including some which inspired this project. 
This work has been supported in part by the U.S. 
Department of Energy under contract No. DOE- 
FG03-97ER40546. 
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